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Application of first and second 
order biodegradation rate kinetic 
model on remediation of crude 


oil polluted soil environment 


Uku Eruni Philip’, Sakwe Adianimovie’ 


ABSTRACT 


This research deals on the application of first and second order biodegradation 
rate kinetic model on remediation of crude oil polluted soil environment. 
Experiments were made on swampy and clay soil environment using moringa 
seed shell and elephant grass which was grounded into powdered form as 
remediants and later compared with NPK fertilizer and yeast as part of the 
bioremediation process. The obtained results were used to determine the 
maximum specific rates, dissociation constant as well as the kinetic values in 
terms of first and second order kinetic. The behaviour of TBC and TPH in clay 
soil were similar to swampy soil. In the model evaluation, the correlation 
coefficient values (R*) obtained from the Michalis-Menten model were higher 
compared to first and second order degradation rate models. Thus, the ranges of 
R? value obtained from the predictive models are 0.9053 to 0.9794 for first order 
degradation rate, 0.7904 to 0.8796 for second order degradation rate and 0.8839 to 
0.9979 for Michalis-Menten Equation. Hence, the Michalis-Menten Equation 
predicted the residual TPH better than the first and second order degradation rate 
models. However, the models are capable of predicting TPHs content in soil 
amended with moringa seed shell in powdered form, yeast and NPK in swampy 
soil (moringa oleifera) and elephant grass in powdered form, yeast and NPK in 


swampy soil (pennisetum purpureum) during bioremediation process. 


Keywords: Bioremediation, degradation, pollution, environment, model, first and 
second order 


1. INTRODUCTION 


The feasibility of bioremediation depends on the location of contaminants. 
Approaches for implementation of bioremediation depend on whether the 
impacted soil to be treated is intact in the environment or it can be excavated for 
treatment in an offsite facility. If on site, the term in-situ remediation suffices and 
if off-site, it is described as ex-situ (Angelucci et al., 2017; Dorado et al., 2015). The 
term bioremediation has been made of two parts: “bios” means life and refers to 


living organisms and “to remediate” means to solve a problem. 
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Bioremediation can also be described as a biological method used for the removal of contaminants from infected or polluted 
environment (Mapscarta, 2020). A point to emphasize here is that bioremediation and biodegradation should not be confused with 
each other. However, in this research, we tend to research on other biological methods possible for bioremediation, hence we 
adopted the use of moringa seed shell, NPK fertilizer and the application of First and Second Order Biodegradation Rate Kinetic 
Model for a better and effect bioremediation process (Nwankwegu et al., 2016; Patel and Kumar, 2016; Pino-Herrera et al., 2017; 
Talha et al., 2018). 


2. METHODOLOGY 


The First Order Biodegradation Rate Kinetic Model 
The first order biodegradation rate kinetic model for prediction of TPH reduction was developed using the principle of mass 


conservation in a batch reactor. A typical batch reactor (Figure 1). 


IPH,.E, Z | 


Z, X 

Figure 1 Batch Reactor 
The principle of mass conservation is stated as 

Inflow of ie of Rate of Rate of 

. mass from ; 
mass into = ne a ; degradation . accumulation (1) 
system eye due to _ of mass 
reaction within system 


The degradation reaction that takes place in the reactor can be represented by equation (2). 


TPH + E+ Z —2—» P+E4Z (2) 
Where; 
TPH =Total hydrocarbon (pollutant) (g) 
E =Bacteria 
Z = Soil (kg) 
p = Products 


k , = Degradation rate constant (unit according to model used) 


From equation (1) we have 


Inflow of mass into system = QO, Cypy(o) (3) 
Outflow of mass from system = QC TPH (4) 
Rate of TPH degradation = —l;p ay (5) 
ae, d(Crpy V) 
Rate of accumulation = — ——————_ (6) 
dt 
Substituting equation (3) through (6) gives 
d(Cyp_V) 
= TPH 
OC ro) = OC ppy — TrpyV + a (7) 
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Since there is no flow of materials in a batch reactor, the inflow of mass into reactor is equivalent to the out flow of mass from 
reactor. 
Thus, 


OC rH) = OC ppy = 0 (8) 


Also, for a batch process, volume of reactor (vessel) is constant hence, the accumulation term is 


dC 
—y Str 0) 
dt dt 


Therefore, equation reduces to 


dC. 
—TppyV =-V— 


dt 


dC. 
Or ~~ lrpy = ae (10) 


Assuming that the degradation of TPH is described by first order kinetics, then, we have 


dc. 
—Trpy =~ = K Copy (11) 


dt 


Where; 

Q > — Inlet volumetric flow rate (kg/day) 

0) = Outlet volumetric flow rate (kg/day) 

C TPH(o) = Initial concentration of Pollutant (TPH) (mg/kg) 

C rpy — instantaneous concentration of Pollutant (TPH) (mg/kg) 
V = Volume of reactor (m3) 

Vrpy = Rate of TPH degradation (mg/kg.day) 


k qd — IPH degradation rate constant (day~!) 


t = Time of TPH degradation (day) 
Integration of equation (11) by the separation of variable method yields the following 


Creat) Z 
dC 
—TPeH —_k, | dt (12) 
Cc. Cren 0 
TPH (0) 
Cc 
In| —“2° | =-k,t (13) 
Creu) 
In Creu) —In Crenco) =—K,t 
In CrpH() =In CrpH(0) = k yt (14) 
Equation (14) can be simply written as 
In Cw = In C, = k yt (15) 
Equation (15) can be compared with the general linear equation of the form 
y=mx+C (16) 
where; 
y=IncC,, 
x=t 
m=k, 
C=InC, 
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A plot of In C_,, against, gives a linear (straight line) graph with slope equivalent to “ k q and intercept equivalent to Inc as 


(t) 
However, to obtain the instantaneous TPH concentration, exponential of both sides of equation (15) is taken to give 
CrpH() = CrpH(0) exp(—k,t) (17) 


Equation (17) is the predictive first order kinetic model for TPH reduction during the treatment process. 


The Second Order Biodegradation Rate Kinetic Model 
The degradation rate of TPH would be further tested with the second order rate kinetics and it is expressed as: 


dC 


2 
Shey maou (18) 
Upon integration of equation (18), we have as follows. 
C; t 
[cvdc =-[k,dt 
C, 0 
1 1 
—-—=k,t (19) 
C, Cy 
Or 
1 1 
— Sth (20) 
C, 0 
Where; 
C 9 and C , are the initial and concentration of TPH at any time, t; while k q is the degradation rate constant for the second order. 
1 
A plot of — against ¢ will give slopeas k, and — as intercept. 
t 0 


Rearranging equation (20) with little manipulations gives the TPH predictive equation using the second order kinetic model as: 


C 
_ TPH (0)t (21) 


TPH ~~ 
ee Ta Crea kot 


3: RESULT AND DISCUSSION 


Evaluation of Treatment Performance in TPH Removal from Soil 

The performances of moringa seed shell and elephant grass in powdered form, yeast and NPK in swampy soil in terms of the 
removal of TPH from polluted swampy soil and clay soils are evaluated (Figures 2, 3, 4). While the effectiveness of the treatments in 
the respective soil was evaluated with respect to time at only 100g weight of treatment. According to Figure 2, the performances of 
the two treatments were compared for each soil with respective to weight at only the 84'" day of the analysis (Figures 3, 4). 

Figure 2 shows the TPH percentage degradation in swampy soil and clay soils obtained at 100g weight of moringa seed shell 
and lephant Grass treatments with time. The percentage of TPH removed from the soils actually increased with time, but the 
removing rate varies for each treatment. The percentage of TPH removed increased from 0% to 97.78%, 97.38%, 98.50% and 97.69% 
for moringa seed shell (SCM), moringa seed shell in clay soil (CM), Elephant Grass in Swampy Soil (GCEG) and Elephant Grass in 
clay soil (CEG) respectively at the 84 day. However, within the first 42 days, the order of the degradation effciency are 
CM>CEG>CSM>CEG, but beyond this period, SCEG performance then slightly edged the others, while the CM which recorded the 


highest performanmce within the first 42 days became the least. 
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Figure 2 Performance of Treatments in Swampy Soil and Clay Soils after 84 Days 
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Figure 3 Comparison of TPH Degradation with CSM and CSEG 


Figure 3 demonstrates the TPH degradation percentage at varying treatment weights in swampy soil. This study revealed that 
TPH degradation percentage increases as treatment weight was increased. However, at the 84™ day of the investigation, elephant 
grass in powdered form, yeast and NPK in swampy soil treatment slightly edged the moringa seed shell in powdered form, yeast 
and NPK in swampy soil treatment in silty clay, though, not significant as depicted. The comparison of the moringa seed shell in 
powdered form, yeast and NPK in swampy soil and elephant grass in powdered form, yeast and NPK in swampy soil treatments in 
swampy soil. However, from the analysis, the TPH degradation percentage in swampy soil at 20g to 100g treatment weight 
increased from 60.79% to 97.78% with moringa seed shell in powdered form, yeast and NPK in swampy soil and 64.95% to 98.50% 
with elephant grass in powdered form, yeast and NPK in swampy soil. 

Figure 4 illustrates the TPH degradation percentage at varying treatment weights in clay soil. Again, it was revealed that TPH 
degradation percentage increased as treatment weight was increased in clay soil, while the elephant grass in powdered form, yeast 
and NPK in swampy soil treatment slightly edged the moringa seed shell in powdered form, yeast and NPK in swampy soil 
treatment. The comparison of the moringa seed shell in powdered form, yeast and NPK in swampy soil and elephant grass in 
powdered form, yeast and NPK in swampy soil treatments in clay soil. Therefore, from the analysis, the TPH degradation 
percentage in clay soil at 20g to 100g weight of treatment increased from 66.89% to 97.38% with moringa seed shell in powdered 
form, Yeast and NPK in swampy soil and 69.74% to 97.69% with elephant grass in powdered form, yeast and NPK in swampy soil. 
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This scenario was observed almost in all the treatment options, although, both treatments proved to be effective for bioremediation 
of Total Petroleum Hydrocarbons content from polluted soil. 
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Figure 4 Comparison of TPH Degradation with CM and CEG 


Degradation Rate Constant for Model Prediction of TPH in Soil 

The degradation of TPH in soil environment has been extensively studied via kinetic models, which has been a useful tool for 
understanding of the reaction dynamics of crude oil during bioremediation (Talha et al., 2018; Tekere, 2019). In literature, the first 
order degradation rate kinetics has been the most applied model for prediction of TPH with time. However, in addition to the first 
order rate model, the second order rate and the Micaelis-Menten models were used in this study. This is to compare the 
predictability of the models and suggest which model best suited the study of TPH degradation in soil amended with moringa seed 
shell in powdered form, yeast and NPK in swampy soil and elephant grass in powdered form, yeast and NPK in swampy soil. To 
effectively utilize these models as expressed the kinetic parameters must be first determined (Ukpaka and Nkakina, 2017). 


First Order Degradation Rate 

The degradation rate constant in the first order rate kinetic model was determined by using equation 4. The plots as in Figure 5, 6, 7 
and 8 illustrates the first order rate kinetics of swampy soil treated with only powder form of moringa seed shell in powdered form, 
yeast and NPK in swampy soil for various reactors conditions of 20g, 40g, 60g, 80g and 100g. The rate values where determined 
based on mass of moringa seed shell in powdered form, yeast and NPK in swampy soil powder form, yeast and NPK introduced in 
each reactor. 


20g: y = -0.0119x + 10.866 
R? = 0.9339 


40g: y = -0.0198x + 10.972 
R2 = 0.9372 80g: y = -0.0358x + 10.97 
R? = 0.9853 


Log (TPH) 


60g: y = -0.0286x + 11.013 
R2 = 0.9729 100g: y= -0.0527x + 11.024 


R? = 0.9527 


40 60 


Time (Days) 


Figure 5 Graph of Log (TPH) in Swampy Soil against Time for Mixture of Moringa Seed Shell in Powdered Form, Yeast and NPK 
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20g: y = -0.0134x + 10.808 
R2 = 0.9522 


xX 


40g: y = -0.0207x + 10.858 


R? = 0.9799 80g: y = -0.0362x + 10.797 


R? = 0.9773 


=x 
a 
= 
{ols} 
fe} 
o_d 


60g: y = -0.0299x + 10.863 


R? = 0.9894 100g: y = -0.0483x + 10.702 


R? = 0.9068 


40 60 
Time (Days) 


Figure 6 Graph of Log (TPH) in Clay Soil against Time for Mixture of Moringa Seed Shell in Powdered Form, Yeast and NPK 


20g: y = -0.0133x + 10.903 
R? = 0.9298 


AOg: y = -0.0221x+ 11.044 
R? = 0.9063 


Log (TPH) 


y60g: = -0.0272x + 11.022 
R? = 0.937 
100g: y = -0.0573x + 11.2 


80g: y = -0.0386x + 11.052 R2 = 0.947 


R? = 0.9549 


40 60 


Time (Days) 


Figure 7 Graph of Log (TPH) in Swampy Soil against Time for Mixture of Elephant Grass in Powdered Form, Yeast and NPK 


Figures 5, 6, 7 and 8 were used to calculate the degradation rate constant, k for the first order rate kinetic model. Thus, from the 
given equations on the plots and comparing with equation 14 the respective predictive model can be obtained for each of treatment 
in the soils. For instance, the respective predictive model determined from the first order rate model for Swampy Soil treated with 
moringa seed shell in powdered form, yeast and NPK in swampy soil (Table 1). For instance, the predictive model for residual 
TPHs in swampy soil treated with moringa seed shell in powdered form, yeast and NPK in swampy soil at the different treatment 
options is presented (Table 1). The degradation rate constant (k) for the respective treatment options (Table 1), which has been 
inserted in the original model. 

From Table 1, it was observed that the degradation rate constant, k increased with increase in the weight of moringa seed shell 
in powdered form, yeast and NPK in swampy soil treatment. The variable x in the regression equation stands for time, while the 
coefficient of x is the TPH degradation constant k. The negative sign in the equation indicates that TPH is reduced with time, and 
R’is the correlation coefficient, which shows the degree of agreement between the experimental value and the predicted 
counterpart. Also, the degradation rate constant, k reported in previous works for first order rate model (Table 2). Hence, the TPH 
degradation constant k obtained in this study are within the reported values in previous works. 
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Figure 8 Graph of Log (TPH) in Clay Soil against Time for Mixture of Elephant Grass in Powdered Form, Yeast and NPK 


2.00 
0.00 


4 
8.00 | agg: y = -0.0242x + 10.974 


R? = 0.9595 
6.00 


60g: y = -0.0334x + 11.026 
4.00 R? = 0.9623 


20g: y = -0.0148x + 10.878 


xX 


R? = 0.9351 


80g: y = -0.0377x + 10.887 


R? = 0.97 


2.00 100g: y = -0.0507x + 10.884 
R? = 0.9456 


0.00 
40 


60 


Time (Days) 


Table 1 Degradation Constant and First Order Rate Model 


Weight (g) | Regression Equation K (day*) | Predictive Model 

20 InC,,, =—0.01191 +10.866 | 0.0119 | Cy, =51067 .73e °°" 
40 In C,p, = —0.0198t + 10.972 | 0.0198 Coy = 51067 736°" 
60 In Crp, = —0.02861 +11.013 | 0.0286 | Crp, =51067.73e 08" 
80 In Cypy, = —0.0358t +10.970 | 0.0358 Crpy = 51067 .73e °°?" 
100 In Cypy = —0.0527¢ +11.024 | 0.0527 | Cypy = 51067 .73e 0?" 


Table 2 Reported TPH Degradation Rate Constant (ka) for different Treatments 
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Treatment ka (day) Reference 
Beans Shell 0.0251 - 

Cassava Peel 0.0288 - 

Cow Dung 0.0498 - 
Groundnut Shell | 0.0260 - 

Melon Shell 0.0257 - 

NPK Fertilizer 0.0228 - 

Pig Dung 0.0266 7 

Cow Dung 0.016 

NPK Fertilizer 0.025 - 

Poultry Manure 0.017 - 

Compost 0.0352 

Spent Mushroom | 0.0366-0.0386 

20g treatment 0.0119-0.0148 | This Work 
40g treatment 0.0198-0.0242 | This Work 
60g treatment 0.0272-0.0334 | This Work 
80g treatment 0.0358-0.0386 | This Work 
100g treatment 0.0483-0.0573 | This Work 
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Second Order Degradation Rate 

The degradation rate constant for the second order rate kinetic model was determined using equation (19) upon obtaining the plots 
(Figure 9, 10, 11, 12), the regression equations were then used to evaluate the rate constant for the different weights investigated in 
the study. 


20 40 60g <80g 100 
g M40g A60g g 8 20g: y = 4E-07x + 2E-05 


yo 
0.001 40g: y = 1E-06x + 7E-06 R°= 0.916 


R? = 0.822 
0.0008 


60g: y = 2E-06x - 1E-05 


100g: y = 1E-05x - 0.0001 


2— 
80g: y = 4E-06x - 3E-05 R? = 0.8761 


0.0004 
R? = 0.8877 x 


0.0002 


0 


Time (Days) 


Figure 9 Reciporical of Lepy in Swampy Soil against Time for Mixture of Moringa Seed Shell in Powdered Form, Yeast and NPK 


Figure 9, 10, 11 and 12 show the plots for calculation of TPH degradation rate constant, k for the second order rate kinetic model, 
which were obtained using with equation (19). Again, from the regression equations on the plots, the second order TPH 
degradation rate constants were evaluated for the different weights of treatments. 


@20g ™40g A60g X80g 100g 20g: y = SE-O7x + 2E-05 
R? = 0.9396 


0.0009 AOg: y = 1E-O06x + 8E-06 

0.0008 | 60g: y = 2E-06x - 8E-06 R? = 0.8478 

0.0007 R? = 0.9063 

0.0006 XK 

0.0005 100g: y = 1E-05x - 7E-05 


80g: y= 4E-06x - 3E-05 R2 = 0.9304 
o ones R? = 0.9234 x 
0.0003 


0.0002 
0.0001 
0 
-0.0001 40 60 
-0.0002 


Time (Days) 


Figure 10 Reciporical of 1 /ppyin Clay Soil against Time for Mixture of Elephant Grass in Powdered Form, Yeast and NPK 
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20g M40g A60g X80g 100g 20g: y = 5E-07x + 2E-05 


0.0014 R? = 0.9068 
4Og: y = 1E-O6x + 2E-06 XK 
0.0012 R? = 0.8004 


0.001 
0.0008 
0.0006 
0.0004 


60g: y = 2E-06x - 3E-06 
R2 = 0.8662 


80g: y = 4E-06x - 4E-05 
R? = 0.8812 


100g: y = 2E-05x - 0.0002 
R2 = 0.8526 


0.0002 


40 60 


Time (Days) 


Figure 11 Reciporical of 1 /ppyin Swampy Soil against Time for Mixture of Elephant Grass in Powdered Form, Yeast and NPK 


@20g M40g A60g <X80g 100g 20g: y = 6E-O7x + 1E-05 
R? = 0.9096 
40g: y = 1E-06x + 1E-07 60g: y = 3E-06x - 2E-05 
80g: y = 4E-06x - 3E-05 * 
0.0006 R? = 0.9294 


0.001 


0.0004 | 100g: y = 1E-05x - 0.0001 
R2 = 0.9284 
0.0002 


Time (Days) 


Figure 12 Reciporical of 1 ppyin Clay Soil against Time for Mixture of Elephant Grass in Powdered Form, Yeast and NPK 
Table 3 TPH Degradation Predictive Model using the Second Order Rate Model 
Weight (g) | Predictive Model 
" — 51067.73 
1+0.0204t 
46 _— 5 1067.73 
Le O.05 1g 
a = 5 1067.73 
1+0.1021¢ 
a be 51067.73 
1+ 0.2043t 
ve bates 51067.73 
1+0.5107¢ 
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The coefficient of time, t in the respective equations (Table 3) is the products of the initial TPH concentration and the second 
order degradation rate constant. The time coefficient increases with increase in the weight of treatment. 


4. CONCLUSION 


The elephant grass treatment slightly edged the moringa seed shell in terms of TPH removal from both soils. This is justified by the 
higher content of total nitrogen and phosphorus in soils amended with elephant grass, which supplied more nutrients to the 
hydrocarbon degrading bacteria. However, both treatments showed high promise for TPH remediation in polluted soils. The 
results obtained from the predictive models, the degree of TPH degradation prediction between the Michalis-Menten rate model 
and the first order degradation rate model cannot be easily distinguished, performed better than the second order kinetic rate. 


Informed consent 


Not applicable. 


Ethical approval 
Not applicable. 


Conflicts of interests 
The authors declare that there are no conflicts of interests. 


Funding 
The study has not received any external funding. 


Data and materials availability 
All data associated with this study are present in the paper. 


KEEFERENGES AND NGIES 


. Angelucci DM, Stazi V, Daugulis AJ, Tomei MC. Treatment of 
Synthetic Tannery Wastewater in a Continuous Two-Phase 
Partitioning Bioreactor: Biodegradation of the Organic 
Fraction and Chromium Separation. J Clean Prod 2017; 152:32 
1-329. 

Dorado AD, Dumont E, Munoz R, Quijano G. A Novel 
Mathematical Approach for the Understanding and 
Optimization of Two-Phase Partitioning Bioreactors Devoted 
to Air Pollution Control. Chem Eng J 2015; 263:239-248. 

. Mapscarta. Places and Locations in Nigeria 2020. https://map 
carta.com 

Nwankwegu AS, Orji MU, Onwosi CO. Studies on Organic 
and In-Organic Biostimulants in Bioremediation of Diesel- 
Contaminated Arable Soil. Chemosphere 2016; 162:148-156. 
Patel BP, Kumar A. Biodegradation of 4-Chlorophenol in an 
Airlift Inner Loop Bioreactor with Mixed Consortium: Effect 
of HRT, Loading Rate and Biogenic Substrate. 3 Biotech 2016; 
6(2):111-123. 


Discovery 59, e71d1266 (2023) 


Pino-Herrera DO, Pechaud Y, Huguenot D, Esposito G, 
Hullebusch ED, Oturan MA. Removal Mechanisms in Aerobic 
Slurry Bioreactors for Remediation of Soils and Sediments 
Polluted with Hydrophobic Organic Compounds: An 
Overview. J Hazard Mater 2017; 339:427-449. 

Talha MA, Goswami M, Giri BS, Sharma A, Rai BN, Singh RS. 
Bioremediation of Congo Red Dye in Immobilized Batch and 
Continuous Packed Bed Bioreactor by Brevibacillus parabrevis 
using Coconut Shell Bio-Char. Bioresour Technol 2018; 252:37- 
43. 

Tekere M. Microbial Bioremediation and Different Bioreactors 
Designs Applied. Biotechnol Bioeng 2019. doi: 10.5772/intech 
open.83661 

Ukpaka CP, Nkakini SO. Crude Oil Remediation using Matlab 
Integrated Agricultural Best Management Practice to 
Improved Soil Nutrients. Pet Petrochemical Eng J 2017; 1(1):1 
01-106. 


11 of 11 


